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Abstract The restricted (equilateral) four-body prob- 
lem consists of three bodies of masses TOi , to-2 and 7713 
(called primaries) lying in a Lagrangian configuration 
of the three-body problem i.e., they remain fixed at the 
apices of an equilateral triangle in a rotating coordinate 
system. A massless fourth body moves under the New- 
tonian gravitation law due to the three primaries; as in 
the restricted three-body problem (R3BP), the fourth 
mass does not affect the motion of the three primaries. 
In this paper we explore symmetric periodic orbits of 
the restricted four-body problem (R4BP) for the case 
of two equal masses where they satisfy approximately 
the Routh's critical value. We will classify them in nine 
families of periodic orbits. We offer an exhaustive study 
of each family and the stability of each of them. 
Keywords: Periodic orbits, four-body problem, sta- 
bility, characteristic curves, asymptotic orbits. 
AMS Classification: 70F15, 70F16 



1 Introduction 

Few bodies problems have been studied for long time 
in celestial mechanics, either as simplified models of 
more complex planetary systems or as benchmark mod- 
els where new mathematical theories can be tested. The 
three-body problem has been source of inspiration and 
study in Celestial Mechanics since Newton and Euler. 
In recent years it has been discovered multiple stellar 
systems such as double stars and triple systems. The re- 
stricted three body problem (R3BP) has demonstrated 
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to be a good model of several systems in our solar sys- 
tem such as the Sun- Jupiter- Asteroid system, and with 
less accuracy the Sun-Earth-Moon system. In analogy 
with the R3BP, in this paper we study a restricted prob- 
lem of four bodies consisting of three primaries moving 
in circular orbits keeping an equilateral triangle config- 
uration and a massless particle moving under the grav- 
itational attraction of the primaries. Here we focus on 
the study of families of periodic orbits. We refer to this 
as the restricted four body problem (R4BP). There ex- 
ist some preliminary studies of this problem in different 
versions, 



Simo 



(1978), 



Leandro 



and Baltagiannis & Papadakis (2011b) studied the equi 



(20061, Pedersen (1944) 



librium points and their stability of this problem. Other 
authors have studied the case where the primaries form 
a collinear configuration. At the time of writing this 
paper we became aware of the paper Baltagiannis & 



Papadakis (2011), where they performed a numerical 



study similar to ours for two cases depending on the 
masses of the primaries: (a) three equal masses and (b) 
two equal masses. It is the second case that our work is 



related to Baltagiannis & Papadakis (2011), although 



we use a slightly different value of the mass parame- 
ter. The reason is the same as the cited authors, of 
having the primaries moving in linearly stable circular 
orbits for a value of the mass parameter less but ap- 
proximately equal to Routh's critical value. By histori- 
cal and theoretical aspects, we use the same letters used 
in the Copenhagen category of the R3BP to denote the 
families of periodic orbits, see Szebehely (1967). The 



families g, /, a, to, r2, 54, ge are similar to those families 
denoted by the same letter in the R3BP, i.e., the fam- 
ily of direct periodic orbits around the mass toi of this 
paper is denoted by the letter g as it was done in the 
Copenhagen category for each family, but the families 
j and j 2 are exclusive of this problem because they do 
not have similar families in the R3BP. Our results con- 
firm and extend four families of periodic orbits found 



2 



in Baltagiannis & Papadakis (2011), such families are 



h, hi h, h- These families correspond respectively 
to the first phases (defined in 3.1) of the families a, j, 
f,goi this paper, however we present 5 new families of 
periodic orbits. We used systematically regularization 
of binary collisions of the infinitesimal with any of the 
primaries by a method similar to Birkhoff 's which per- 
mit us to continue some of the families beyond double 
collisions. In this way we can show that such continued 
families end up in a homoclinic connection. This last 
phenomenon can be dynamically explained by the so 
called blue sky catastrophe and a rigorous justification 
will appear elsewhere. 

We recall that Routh's criterion for linear stability 
states that 

777.17712 + 77727773 + 77137711 1 

rrii + m 2 + m 3 27 

When the three masses are such that 777,2 = m 3 := /i 
and mi + TU2 + 7773 = 1, the inequality is satisfied in 
the interval /j € [0, 0.019063652805978857 . . .), so in our 
case we take the masses equal to m 2 = m s = 0.0190636 
and mi = 0.9809364. 



2 Equations of Motion 

Consider three point masses, called primaries, moving 
in circular periodic orbits around their center of mass 
under their mutual Newtonian gravitational attraction, 
forming an equilateral triangle configuration. A third 
massless particle moving in the same plane is acted 
upon the attraction of the primaries. The equations 
of motion of the massless particle referred to a synodic 
frame with the same origin, where the primaries remain 
fixed, are: 



x" — 2ny' — n 2 x 



= -* 2 £'< 



(x - Xi) 



H „3 



(y - Vi) 



y" + 2jix'~n 2 y = -fc 2 ^m /* (1) 
where k 2 is the gravitational constant, n is the mean 



motion, 



(x — x^ + (y — yi) is the distance of 



the massless particle to the primaries, Xj, y~i are the 
vertices of equilateral triangle formed by the primaries, 
and (') denotes derivative with respect to time t* . We 
choose the orientation of the triangle of masses such 
that mi lies along the positive x-axis and 7772, 7773 are 
located symmetrically with respect to the same axis, 
see Figure [T] 

The equations of motion can be recast in dimension- 
less form as follows: Let L denote the length of triangle 



formed by the primaries, x — x/L, y — y/L, Xi = Xi/L, 
Hi — Vi/ L, for i — 1,2, 3; M = mi + 7772 + 7773 the total 
mass, and t — nt* . Then the equations |l]) become 



x — 2y — x 
y + 2x-y 



•i=i 
3 



(X - Xi) 
' i 

(y - yd 



(2) 



i=l 



where we have used Kepler's third law: k 2 M = n 2 L 3 , (') 
represents derivatives with respect to the dimensionless 
time t and rf = (x — x^ 2 + (y — yi) 2 . 

The system ^ will be defined if we know the vertices 
of triangle for each value of the masses. In this paper 
we suppose /i := 7772 = 7773 then n\ = \ — 2/7. It is 
not difficult to show that the vertices of triangle are 
given by X\ = \/3/i, yi = 0, x 2 = - 3 ^T 2 ^ , y 2 = -\, 
X3 = — ^ /3 ~( 1 ~ 2 f J -) ; y 3 — 1 The system (12 h can be written 
succinctly as 



x — 2y = Q x 
y + 2x = fly 

where 



(3) 
(4) 



i 3 

at* +y ) + 

i— 1 



is the effective potential function. 
There are three limiting cases: 

1. If fi = 0, we obtain the rotating Kepler's problem, 
with mi = 1 at the origin of coordinates. 

2. If /7 = 1/2, we obtain the circular restricted three 
body problem, with two equal masses 7712 = 7773 = 
1/2. 




Fig. 1 The restricted four-body problem in a synodic sys- 
tem 
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3. If /j = 1/3, we obtain the symmetric case with three 
masses equal to 1/3. 

It will be useful to write the system ^ using com- 
plex notation. Let z = x + iy, then 



z + 2iz = 2 



with 



dz 



1 



(5) 



n(z,z,fj,) = ^\z\ 2 + U(z,z,fi) 
where the gravitational potential is 

3 



U(z,z,(i) = ^2 



Hi 



i=l ' 

\z - zA, i 



and Ti = \z — Zi\, i = 1,2,3 are the distances to the 
primaries. System ^ has the Jacobian first integral 



2{l(z,z,(x)-\zf 



C. 



If we define P = p x + ip y , the conjugate momenta 
of z, then system (J3j) can be recast as a Hamiltonian 
system with Hamiltonian 



H 



-\P\ 2 +Im(zP)^U(z 1 z, f ,) 

\{pl+pl) + {Wx-xp y )-U{x,y,n). (6) 



The relationship with the Jacobian integral is H 
—C/2. The phase space of ^ is defined as 



A = {(z,P)e 



X «J z 



j^z h i= 1,2,3}, 



with collisions occurring at z = Zi, i = 1, 2, 3. 

There exist five equilibrium points for all values of 
the masses of the primaries in the R3BP; in the R4BP, 
the number of equilibrium points depends on the par- 
ticular values of the masses. For the value of the 
mass parameter we are using throughout this paper of 
fi = 0.0190636, the Hill's regions are shown in Figure[2j 
For large values of C the Hill's regions consist of small 
disks around the primaries together with an unbounded 
component having as boundary a closed curve around 
the primaries. As the Jacobian constant decreases, the 
evolution of the Hill's region is shown in Figure [2] The 
smaller value of C is just above the critical value where 
the Hill's region is the whole plain minus the positions 
of the primaries. 

A complete discussion of th e equilibrium points 



and bifurcations can be found in |Delgado fc Alvarez- 



Ramirez ( 


2003 


), 


Meyer 


(1987 


, |Leandro| (|2006| 


giannis & Papadakis 


(2011b 


)■ 


Simo 


(1978 


). In ( 
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Fig. 2 Hill's regions for a large value of the Jacobian con- 
stant (top-left). Hill's regions for critical values of the Jaco- 
bian constant (top-right and second row). The Hill's regions 
of the last row correspond to a slightly larger value than the 
critical one for illustrations purposes. 

equilibrium points. We use the notation shown in Fig- 
ure [3] for the eight critical points. All of them are un- 
stable except the non-collinear L7 and L%. 



3 Symmetric periodic orbits 

In what follows we consider symmetric periodic orbits, 
i,e., periodic orbits symmetric with respect to the syn- 
odical x-axis, namely orbits which are invariant under 
the symmetry (z,t) — > (z,—t). Thus a symmetric pe- 
riodic orbit is defined by two successive perpendicular 
crossings with the cc-axis. We use as staring point of 
the continuation either nearly Keplerian circular orbits 
for large values of the Jacobian constant or small Li- 
apunov's orbits emerging form L\. It was shown in 
Leandro (2006) that the equilibrium point L\ is unsta- 



ticular problem we have 2 collinear and 6 non-collinear 



ble with eigenvalues ±iw, ±A where uj and A are real 
numbers, so we can use the Liapunov's Center Theorem 
to find periodic orbits around the equilibrium point L\ 
and continue them. 

In the R3BP the equilibrium points L4 and L5 are 
limits of families of periodic orbits of the family g of 
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Fig. 3 The eight equilibrium points (green) for two equal 
masses. 

the Copenhagen category. This phenomenon is known 
as the "blue sky catastrophe" termination principle. 



parameter [ai,a^] for which the family is followed and 
we want to extend it, then 

1. The family remains in itself, i.e., the characteristic 
curve is a closed curve. We call this family a closed 
family. 

2. For a < a± and a > a.2 the family has a natural 
termination for which one of the following amounts 
grow without limit 

• The dimension D of the orbit, defined as the max- 
imum distance to the origin. 

• The parameter a. 

• The period of the orbit. 

The second case is called an open family. 
Note that the principle of termination of a family of 
periodic orbits mentioned in [T] is a particular case of 
the above definition because the period of the orbit in- 
creases without bound. 



In the papers by Buffoni (1999), Meyer & McSwiggen 3.1 The search for periodic orbits 



(2002) and Simo (1978) a complete discussion of this 



principle in the R3BP is given. We state the theorem 



behind this phenomena as stated in Henrard (1973): 
Theorem 1 Let us consider a non- degenerate homo- 
clinic orbit to an equilibrium with eigenvalues ±a ± ito 
with a and u) reals and strictly positive, of a real ana- 
lytic Hamiltonian system with two degrees of freedom. 
Close to this orbit there exists an analytical family of 
periodic orbits with the following properties 

1. The family can be parametrized by a parameter e in 
the interval < e < 5 with S small enough. Let us 
write it as x(e, t) = x(e, t + T(e)). 

2. For e = we have the homoclinic orbit. 

3. The period T(e) increases without bound when e goes 
to zero. 

4. The characteristic exponents of the family change 
from the stable type to unstable type an vice-versa 
infinitely many times as e goes to zero. 

In a future work we will discuss the theoretical aspects 
of this termination principle in the R4BP. 

In the following section we will classify a large num- 
ber of periodic orbits in sets called families. We will 
take the Henon's definition for a family of periodic or- 



The periodic orbits were calculated in double preci- 
sion with a multi-step Adams-Bashforth integrator of 
variable order for more accuracy. New transformations 
were needed to regularize different kind of collisions ap- 
pearing in the families of periodic orbits of this prob- 
lem. The families have been identified by letters as in 
the Copenhagen category with or without subscripts, 
the subscripts meaning the number of loops of the or- 
bit around the primary under consideration. We use 
the classic (x, C) plane of characteristic curves to rep- 
resent the families of periodic orbits, in addition we use 
the (a, C) plane to show the evolution of the stability 
of the families, here a denotes the stability index (see 



Henon (1965b) for details), we have stability in the lin- 
ear sense when \a\ < 1 and instability in other case. 

The families have been separated in phases as in the 
R3BP (see Szebehely (1967)), the colors in the charac- 



teristic curves of the nine families represent the differ- 
ent phases (and orbits near to collision) of each family, 
representative orbits are shown to illustrate each phase. 



Some orbits shown by Broucke ( 1968 1 and Baltagiannis 



& Papadakis (2011) can be compared with ours. 



bits to make such classification (see Henon (1997)): 
Definition 1 A family of periodic orbits is a set of 
symmetric periodic orbits for which the initial parame- 
ter Xq = x(0) and the period in family can be considered 
as two continuous functions of one single parameter a. 

In general we will consider the Stromgren's termina- 
tion principle to decide when a family ends. 
Definition 2 Suppose we have obtained a finite sec- 
tion of a family of periodic orbits in an interval of the 



4 Classification of families of periodic orbits 

4.1 The family g of direct periodic orbits around m\ 

The first phase of this family stars with infinitesimal di- 
rect circular periodic orbits around mi, the size of the 
orbits increases as the value of Jacobi constant C de- 
creases until a collision orbit is reached, the first phase 
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Fig. 4 The characteristic curves of nine families of periodic 
orbits for the restricted four-body problem with two equal 
masses. 

Familiy g 




Fig. 5 Characteristic curve of family g. 



of this family was established and studied in Baltagian- 
nis & Papadakis (2011) for a little smaller value of \i. 



Second phase starts with the retrograde orbit following 
the collision orbit which forms two loops as shown in 
Figure 17 The inside loops now increase their size 



and the outside loops shrink as the value of C deceases, 
both sets of loops become indistinguishable at the fold 
point when C s» 0.1797. After this point the third phase 
starts, now the role of the loops is interchanged i,e., the 



inside loops shrinks and the outside loops expand as the 
value of C increases. 

Following the evolution of this phase we found that 
an orbit of collision with m\ appears and the inside 
loops disappear, this is the beginning of fourth phase 
where the middle part of the orbits increases its size as 
the Jacobi constant increases, the termination of this 
phase (and of the whole family) are asymptotic orbits 
to Z/2 (£3 in Baltagiannis & Papadakis (2011)). More 
precisely, the value of C oscillates in a small neighbour- 
hood around the value of the Jacobi constant of the 
equilibrium point L2 and the period tends to infinity 
as is predicted in theorem [l] similar to the Copenhagen 
category of R3BP. 

4.2 The family / of retrograde orbits around mi 

Family f 
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Fig. 6 Characteristic curve of family /, zoomed area indi- 
cates the end of family at the equilibrium point L2. 

The first phase of this family starts with infinitesi- 
mal retrograde circular periodic orbits around mi, as in 
family g the size of the orbits increase as the value of C 
decreases monotonically until a fold point is reached, 
this happens at C ~ —0.6379. This is the beginning of 
second phase, the periodic orbits still continue increas- 
ing their size but now these orbits tend to collision with 
the primaries mi and 013 , however this collision is never 
reached because the periodic orbits become asymptotic 
to L 2 , as in family g. See Figures |6| [l~8| 

4.3 The family a of retrograde orbits around L\ 

The beginning of this family is provided by the Lia- 
punov's center theorem, therefore in this family we con- 
tinue retrograde periodic orbits around the equilibrium 




Fig. 7 Characteristic curve of family a. 



Fig. 8 Characteristic curve of the closed family 



point L\ (Li in Baltagiannis & Papadakis (2011)) for 
values of the Jacobi constant less than C\. As the Ja- 
cobi constant decreases monotonically, the periodic or- 
bits increase their size until a collision orbit with wi\ is 
reached; this is the end of first phase. While the value of 
C continues decreasing, a second loop appears in the or- 
bits, this loop increases its size along this second phase 
until a new fold point is reached when C « —0.5846. 
At this point the inner and outer loops become indis- 
tinguishable as in the previous families, after this fold 
point, both loops are interchanged and the new inside 
loop shrinks as the value of C increases until a collision 
orbit with ra\ finishes the third phase, see Figures [7J 

m 

4.4 The family g 4 

This family is the first example in this paper of a closed 
family of periodic orbits, see Figures [HJ [20j The first 
phase of this family is composed of periodic orbits form- 
ing four loops around mi, following the evolution of 
first phase we found that as C decreases the four loops 
shrink until they become indistinguishable i,e; a fold 
point is reached at C « 1.3381, the second phase starts 
when these loops separate each other, as C increases 
the periodic orbits tend to collision with mi i,e; the 
intersections between the loops tend to mi. After this 
collision, the third phase starts. The periodic orbits 
change multiplicity because four inside loops around 
mi appear together with four outside loops. As the 
value of C increases, the inside loops increase their size 
while the outside loops shrink at same time until they 
disappear. This is the end of third phase. 



Fourth phase begins when the loops around mi in- 
crease their size as C continues increasing monoton- 
ically. These loops become indistinguishable at the 
value C s» 3.7581 and a new fold point is reached. Af- 
ter this fold point the fifth phase begins, as expected 
inside and outside loops interchange and the peak of 
the outside loops become non-smooth i,e; no more or- 
thogonal intersections with rr-axis exist, this is the end 
of fifth phase. As C decreases the loops of the orbits 
shrink to collision with mi and the sixth phase ends. 
We observe that the first phase of this family starts 
after this collision, therefore the family is closed. 

4.5 The family g$ 

This family is another example of a closed family of 
periodic orbits. The behavior of its orbits is more com- 
plicated than that of previous families see Figure [9j 



21 



We have named g§ this family because there exist 
a section of this family where the orbits form 6 loops 
around the primary mi but following the evolution of 
this family, we find that the orbits show a complicated 
behavior, therefore is not clear at all how to classify in 
phases this family, for simplicity we use the term "first" 
phase to the outer section of the characteristic curve be- 
tween the return points in Figure [9j and second phase 
to the inner section between the return points. In each 



case we show representative orbits in Figure 21 



4.6 The family m of retrograde orbits around mi, ni2 
and m3 

Family m consists of retrograde periodic orbits around 
the three primaries but these orbits do not surround 
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4.7 The family j of retrograde periodic orbits around 
771.2 and 7773 



Fig. 9 Characteristic curve of the closed family g@. 

Familiy m 




Fig. 10 Characteristic curve of the closed family m. 

the primaries simultaneously i,e., they form three loops, 
each loops surrounds one primary as can be seen in the 
Figure (22). The characteristic curve of the family 



is shown in Figure 10 The three loops increase and 
decrease their size while the family is followed; however 
collision with the primaries never is reached although 
the orbits are close to collision when the loops decrease 
their size. This behavior is cyclic because the family is 
closed. 



Familiy j 
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Fig. 11 Characteristic curve of the family j. 

The first phase of this family by one side tends to col- 
lision with the primaries 7712 and 7773, such collisions are 
reached for great values (negatives) of the Jacobi con- 
stant, see Figures [TT] [23} As C increases, the orbits 
increase their size around 7772 and 7773, the orthogonal 
intersection to the right of both primaries 7772 and 7773 
tends to the primary mj and a collision orbit appears, 
this is the end of the first phase. After this collision an 
inner loop appears as expected and the orbits become 
direct around m\, this is the second phase, the inside 
and the outside loops increase their size and the orbits 
tend to collision with both primaries 7772 and 7773 until 
this collision happens. The third phase starts after this 
collision, two loops around 777.2 and 7773 respectively ap- 
pear in the orbits and as C decreases these loops these 
loops increase around the primaries, the behavior of the 
orbits complicate while the family is followed. 

Finally a new collision orbit with mi appears. We 
have decided to terminate the family at this point, the 
complicated behavior of the orbits and the long time 
of integration of the regularized equations forced us to 
stop the continuation at this point. 

4.8 The family ri of asymptotic orbits to L2 

This family of periodic orbits has been named ?"2 in 
analogy with the family r of the Copenhagen category, 
the subscript indicates that the family is asymptotic to 
1/2 see Figures [T2j [24] The first phase of this family is 
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Family r. 



° -2 




Familiy j2 



Fig. 12 Characteristic curve of family r2, zoomed areas 
indicate the end of family at the equilibrium point L2 and 
a return point of the family. 



composed by asymptotic periodic orbits to the equilib- 
rium point L2, these orbits form two loops surrounding 
the primaries and 777,3 but while the periodic or- 
bits go away from Li these loops shrink, therefore we 
have orbits close to collision now, however such colli- 
sion never is reached. When the value of C begins to 
decrease monotonically these loops increase its size and 
therefore the period of the orbits increase, it is interest- 
ing to note that the orbits become symmetric respect 



to y-axis (see Figure 24 ) however this symmetry disap- 



pears as C continues decreasing, a collision orbits with 
mi terminates the first phase. 

As expected, a new loop around mi appears in the 
orbits, here the second phase starts, following the evolu- 
tion of the family we can see that the resulting loops of 
the collision increases its size as C decreases monotoni- 
cally and the orbits tend to collision with the primaries 
7712 and 7773 . At C ps —0.7726 we find a return point. 
We must emphasize that such collision is not reached 
although the mentioned loop continues increasing its 
size and therefore the orbits increase their period. This 
behavior continues as C increases monotonically. As 
in family j the long period of the orbits (long time of 
integration of the equations is needed) and the high in- 
stability of the orbits did not allow us to establish the 
end of this phase and therefore of the whole family. 



0.5 




Fig. 13 Characteristic curve of the family j'2, note it is 
very similar to the one of the family j. 



4.9 The 

7B2 and 7773 



family j 2 of retrograde periodic orbits around 
md 7773 



This last family named J2 contains retrograde periodic 
orbits around 7772 and 7773 as in family j but this time 
the orbits have an extra loop around 7712 and 7773 and 
this loop surrounds the primary mi too (see Figures 13 
25). The evolution of this family is very similar to the 
evolution of family j, in fact, its characteristic curve 
has the same form that the one of family j. We show 
representative orbits of each phase in Figure [25} 



5 Stability of families and critical points 



5.1 Critical points 



As can be seen in Henon (19971, when we follow fam 



ilies of periodic orbits it can be found some especial 
points in such families called critical points such as fold 
points and branching points, the last point is where two 
characteristic curves intersect, in 5.2 we will see that 
these points are in relation with stability changes, this 
facts do not ocurr by casuality, such phenomena was 



studied by example in Henon (1965b). We must say 



that some families of periodic orbits were found through 
these points, in Figure [l4]we show the branching points 
found in this study. 

1. Family g has 3 branching points; first one P\ at C pa 
3.7581 with family #4, second one P2 at C pa 2.1662 
with family g^ again, third one P3 at C pa 0.1797 
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Fig. 14 Branching points between the families of periodic 
orbits 

with family /. This last branching point happens at 
a fold point. 

2. Family / has 5 branching points, first one P4 at C ~ 
3.6364 with family 34, second one P5 at C « 2.8481 
with family <?, third one Pq at C ~ 1.3381 with fam- 
ily (74, fourth one P3 with family g already men- 
tioned, fifth one P7 at C ~ —0.5846 with family a 
(fold point of family a) 

3. Families j and ji intersect at two points, first one P$ 
at C w —1.5015 and second one Pg at C « —1.073. 

5.2 Stability of families 

In family g we see that the first and second phase con- 
tain stable periodic orbits, at the end of second phase 
in the critical (fold) point we have that \a\ = 1 as 
expected, after this critical point all orbits become un- 
stable, at the end of this family we observe strong os- 
cillations of the sign of a and therefore between the 
stable and unstable areas as is predicted by the "blue 
sky catastrophe" termination. In family / we have that 
the first phase of this family posses stable orbits as is 
shown in the Figure |15[ at the fold point we have 
I a I = 1 again, after this point all orbits are unstable un- 
til the termination of family is reached and as in family 
g strong oscillations between stable and unstable orbits 
are observed. 

In family a almost all orbits are unstable but we 
have 3 small regions where the orbits are stable, one of 
them is in the first phase, the second one is at the end 
of second phase; in the fold point we have of course 
\a\ = 1. After this point the stability index starts to 
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Fig. 15 Stability curve of family g (up left), stability curve 
of family / (up right), stability curve of family a (down left), 
stability curve of family g& (down right). 
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Fig. 16 Stability curve of family j (up left), stability curve 
of family 32 (up right), part of the end of the stability curve 
of family r-2, (center). 

decrease and therefore we find a stability region at the 
beginning of third phase. The family 34 is a closed 
family, a "half of its orbits are stable and a "half are 
unstable, this affirmation can be seen in Figure [l5j at 
the two fold points we have \a\ = 1 and between these 
two points the stability changes occur. 

In family g e almost all orbits are unstable, only a 
very little region around the two fold points has stable 
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orbits. Family m presents a similar behavior as in fam- 
ily gg. Orbits in family r 2 are strongly unstable except 
at the termination of family where oscillatory behavior 
between areas of stability and instability is observed. 
In the family j we found a stability region between the 
values C ~ —1.27 and C ps 2.59, see Figure 16 in this 
interval we have 3 critical points on the characteristic 
curve, and this can be seen on the stability curve where 
we have values for which \a\ = 1. For the family ji we 
have its stability curve behaves similar to the one of the 
family j, see Figure 



1G 



6 Examples of periodic orbits 



In Figures [T7j[25]we show representative orbits of each 
phase of the nine found families of periodic orbits. Each 
row in the following figures represent a phase of the 
family, the first row represent the first phase of the 
family, second row represent the second phase etc. 
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Fig. 17 Phases of family g. 



7 Conclusions and remarks 



In this work we extend the previous study of [Baltagian- 
nis &: Papadakis| ( |2011[ ) and we include 5 new families of 
periodic orbits. We performed the numerical continua- 
tion beyond different collisions with the primaries and 
show how some of them end up in a homoclinic orbit 
as predicted by the "blue sky catastrophe" termination 
principle. In the former section we show representa- 
tive orbits for each phase of the families, the orbits are 
mainly retrograde but there are families consisting of 
direct periodic orbits. We have also studied the stabil- 
ity of each family of periodic orbits and we show the 
characteristic curve of seven families where a consider- 
able number of stable periodic orbits were found. Our 
main results can be summarized as follows: 

1. The families g, f , a, g^, j, j2, and T2 present a large 
number of stable periodic orbits, all the orbits in the 
families ge and m are unstable except in very small 
regions. 

2. The families g, f and ri are asymptotic to the equi- 
librium point L2, i,e; they end according to the "Blue 
Sky Catastrophe" termination principle. 

3. The families (74, ge and m are closed families. 

4. The families g, f, a, g<±, j, and j'2 present Branching 
points between their characteristic curves. 

5. Almost all the families consist of retrograde periodic 
orbits, except families g, (74 where have direct ones. 

Acknowledgements Author Burgos-Garci'a has been 
supported by a CONACYT fellowship of doctoral stud- 
ies. 




Fig. 18 Phases of family f. 





Fig. 19 Phases of family a. 
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Fig. 22 Phases of family m. 



Fig. 25 Phases of family jr 2 . 
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